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Abstract 



Three- and four-neutron systems are studied within the framework of the 
hyperspherical approach with a local 5-wave nn-potential. Possible bound 
and resonant states of these systems are sought as zeros of three- and four- 
body Jost functions in the complex momentum plane. It is found that ze- 
ros closest to the origin correspond to sub-threshold (nnn) ^ and (nnnn) 
0"'" resonant states. The positions of these zeros turned out to be sensitive 
to the choice of the nn-potential. For the Malfliet-Tjon potential they are 
Eo{^n) = -4.9 - i6.9 (MeV) and Eoi'^n) = -2.6 - i9.0 (MeV). Movement of 
the zeros with an artificial increase of the potential strength also shows an 
extreme sensitivity to the choice of potential. Thus, to generate ^n and ^n 
bound states, the Yukawa potential needs to be multiplied by 2.67 and 2.32 
respectively, while for the Malfliet-Tjon potential the required multiplicative 
factors are 4.04 and 3.59. 



I. INTRODUCTION 

The question of the existence of pure neutron nuclei is a longstanding problem, almost 
as old as nuclear physics itself. Many experiments have been performed in the past in order 
to find di-, tri-, or tetra-neutron bound states. On the theoretical front much effort has 
also been devoted to understand the physics of few-neutron systems. 



Early experimental and theoretical investigations were controversial and their findings 
varied from non-existence to 'discovery' and 'prediction' of such states. Nowadays it is gen- 
erally accepted that no few-neutron nuclei exist. However, the interest in pure neutron 
systems did not wane. Since (np), (nnp), (npp), and {nnpp) systems are bound, one could 
expect that due to isotopic invariance a small increase in the neutron-neutron attraction 
could generate bound states in their isotopic partners (nn), (nnn), and (nnnn). Indeed, 
this is the case for the two-neutron system which requires only an additional 67keV at- 
traction for the formation of a bound di-neutron. For three and more neutrons, however, 
the situation is more complicated due to the Pauli principle. This becomes clear when one 



*Permanent address: Joint Institute for Nuclear Research, Dubna, 141980, Russia 



1 



considers a multi-neutron system in even the simplest model, namely, that of neutrons mov- 
ing independently in a harmonic oscillator well. Since only two neutrons can occupy the 
ground S'-wave state of this well, the third, fourth, and other neutrons must be in excited 
states and such a multi-neutron system can form quasi-stationary states instead of bound 
ones. This encouraged research by both experimentalists and theorists to locate resonances 
in few-neutron systems 

In spite of the abundance of experimental works that exploit different reactions and 
methods, the situation concerning few-neutron resonances remains unclear. In many cases 
certain irregularities in cross-sections are observed, which can be attributed to such reso- 
nances, but their positions and widths have not been firmly established. As for theoretical 
studies of few-body resonances, they are complicated and cumbersome to the extend that 
these states have not been properly investigated. Since bound states can be located more 
easily, it is customary for theorists to increase the nn-potential until a bound multi-neutron 
system is supported. It was found that the interaction needs to be increased by a factor 
of between 2 and 4, depending on the number of neutrons, the type of potential, and the 
method used. 

There are many different definitions of quantum resonances as well as many methods of 
locating them Though all methods are physically justified and able to produce reason- 

able results, only few of them can be considered as rigorous. These are the methods based 
on the correspondence between resonances and singularities (poles) of the 5'-matrix in the 
unphysical sheet of the complex energy plane. Apart from being mathematically rigorous, 
such approaches enable one to obtain the position and width of a resonance simultaneously 
with the same accuracy, simply as coordinates of a pole on the energy plane. And again, 
there are several different methods for locating the S'-matrix poles (a survey of them can be 
found in Refs. However due to the existence of redundant poles in the S-matrix it is 

preferable to search for zeros of the Jost function instead. 

For two-body systems the Jost function is well-defined and an exact method of calculat- 
ing it has been developed However, for many-body systems, no such method exists 
and in fact no rigorous definition of the Jost function has been given in a general form yet. 
Only for a special class of many-body systems has a definition and method of calculating 
the Jost function been proposed [0, namely, for systems which cannot form clusters of any 
type and thus only one allowed channel exists. Some authors call them 'democratic systems' 
0. Wave functions which describe such systems behave asymptotically as diverging hyper- 
radial spherical waves . Clearly a multineutron system is one such system and it may be 
investigated via the hyperspherical approach. 

In this work we consider three- and four-neutron systems in the minimal approxima- 
tion, L = Lmini of the hyperspherical harmonic approach with local S'-wave nra-potentials. 
Possible bound and resonant states are sought as zeros of three- and four-body Jost func- 
tions in the complex plane of the momentum conjugate to the hyperradius. These zeros are 
related to the total energy in a straightforward way. The many-body Jost function for the 
hyperradial equation is defined similarly to the two-body case and is calculated by a method 
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which combines the variable constant and complex rotation methods 0. The zeros found 
represent the following sub-threshold resonant states: ^n{^ ), '^n{0'^), ^n(l+), and ^n(2"'"). 

Our paper is organized as follows: In Sec. II we discuss the many-body Jost function 
and in Sec. Ill we give the required matrix elements of the potential. In Sec. IV the complex 
coordinate rotation method is used to achieve an analytic continuation of the Jost function 
into the lower half of the momentum plane, while in Sec. V we describe the physical input 
used in the calculations and present our results together with some discussions. 



II. MANY BODY JOST FUNCTION 

Let {r} and {p} be the complete sets of Jacobi vectors defining a multi-neutron con- 
figuration in coordinate and momentum space respectively. The wave function \E'"(r,p) of 
an A^-neutron system can be expanded in terms of hyperspherical harmonics Y^1^{Q) with 
unknown coefficient-functions (r, p) 

oo 

^^{r,p)=r'~'^/' Y: r[2](^)^[L](^,P), (1) 

where {a} is the set of all conserved quantum numbers (total angular momentum, isospin, 
and parity) and the multi-index [L] represents the 'grand orbital' L and other nonconserved 
quantum numbers; r and Q, {r,Q} = {r}, are the hyperradius and hyperangles. For two 
neutrons, N = 2, the expansion is the usual partial-wave decomposition where [L] = [i, m] 
is the angular momentum and its third component, and imin = 0. When N > 2, L starts 
from a generally nonzero value L^in defined by the symmetry properties of the state under 
consideration [|TTI|. 



Substituting the expansion (|T]) into the Schrodinger equation, one arrives at the following 
infinite system of coupled hyperradial equations 

92 + p2 _ ^ ^)/^2j ^a^^^^^ ^ ^ 1^[1][L'](0 U^L']ir, p) • (2) 

IL'] 

Here [L] = [Lmin], [Lmin + 1], ■ ■ ■ and p may be called 'hypermomentum' since it is related to 
the total energy E and the neutron mass m in the same manner as in the two-body case, 
= 2mE (we use units such that ^ = 1); A is an analog of the orbital angular momentum, 

A^L + ^(7V-2), (3) 

which assumes half-integer values for uneven particle number A^. The potential matrix W 
has elements 

= 2m(y[2]l E y^A^)\Y[l']) , (4) 
where the integrals are over the hyperangles and Vij are the two-body potentials. 
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The required boundary conditions for the differential equations (0) are, firstly, the solu- 
tions must be regular at the origin, 

«f^](r,p)^0, r^O, (5) 

and, secondly, at infinity must be of some physically motivated form 

uU^,P)-^U^L]ir,p) , r^oo. (6) 

In the general case the boundary function [/['^] (r, p) consists of terms describing each open 
channel. Since a multi-neutron system is a 'democratic' one, we shall construct t/j*}^] as an 
one-channel boundary function. 

Instead of matching u'^i^-^ir, p) to U^i^^{r, p), we can secure the one-channel boundary con- 
dition automatically, and in an exact form, by the following procedure. Let us consider a 
general regular solution of (H), defined only by the condition (H) and without any restric- 
tions at large r. Certainly many such solutions exist, each having different behaviour at 
large distances. By choosing only those which are linearly independent, we have a regular 
basis. Any other solution which is regular and has a specific behaviour at large r must be 
a linear combination of this basis set. Thus, we can find first the regular basis and then 
construct the required physical solution. 

Like any other basis, the regular basis may be chosen in an infinite number of ways. 
Exploiting this freedom, we can choose the basis regular solutions in a way which is most 
suitable for the subsequent construction of a physical solution. Considering the system (H) 
(which is always truncated at some [Lmax]) as a matrix equation, we see that it has as many 
independent regular column-solutions as the column dimension (number of equations in the 
system). We may combine all these linear independent columns in a square matrix 
of the form (for simplicity we drop the superscript a ) 

0[L][L'](r,p) = ^ {h[^\pr)Fl^l,^ir,p) + h[-\pr)Fl-l,^{r,p)} , (7) 
where h^^^ are the Riccati-Hankel functions [Jl2[, and F!"^hn{r,p) are new unknown func- 



tions. Since instead of one unknown matrix we have introduced two matrices F|^^j]^>■^, 

we require an additional constraint. The most convenient is the Lagrange condition 

ipr)drFl^l,^ (r, p) + h[-^ (pr)a.F^,] (r, p) = , (8) 

which is standard in the variable-constant method for solving differential equations ||T3| . 



Substituting the ansatz (0) into Eq. (H) and using the condition (|), we derive the 
following coupled differential matrix equations of first order: 

drFl^]lL']ir,p) = E W^[L][L"](r) {4"^(pr)FWj^,,(r,p) + ipr)Fl-\,,^{r,p)} , 

(9) 

fLUr,p) = E ^[L][L"](r) {4t^(pr)4+;,,,(r,p) + h[-) {pr-)Fl[}^^^,^{r^p)] . 



drF, 
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These equations must be supplemented with appropriate boundary conditions at r = 0. 
In Ref. [14 1 it was shown that for an arbitrary A'"~body system the fundamental system of 
regular solutions of Eq. (0) vanishes near r = in such a way that 



(p[L]lL']{r,p) _ , . 



Thus, we can define the regular basis by the following boundary condition 

0[L][L'](r,p) _ , . 

where jx is the Riccati-Bessel function [|12|. This is in accordance with (|1^) and is a natural 
generalization of the corresponding boundary condition of the two-body problem. 

Since, by definition, is regular at r = 0, the behaviour of the functions F^]j\il'] 

and F^L][L'] of Eq. near the origin is such that the singularities of h^x\p^) ^'^d \pr) 
compensate each other. This can be achieved if and F^^^ are identical as r — 0, i.e., 

for then 

<P[L][L']{r,p) ^r^^jxipr) FiL]iL,^{r,p) , (12) 



^[L][L'] ■ (13) 



so that the boundary conditions, Eq. (|TT]), become 

jx{pr)F^^lj^,^{r,p) 



lim 



Jy{pr) 



However, since r = is a singular point, for practical calculations one needs to solve the 
system (^ analytically on a small interval (0, 6] and then impose the boundary conditions 
at r = 6. Such an analytical solution can be easily found by choosing 6 to be small enough 
so that for r G (0, 6] we may write 

drFl^l']{r,p) = ±^h^Pipr) W^[L][L'](r) jx'ipr) • 
tp 

For small r the Riccati-Neumann function nx is dominant in h^x^^ = jx ± inx and thus 

9rFl^lL,]{r,p) ^ --nx{pr) W[L][L']{r) jx'{pr) . 

Upon integrating this (approximate) equation we find 

1 r 

^[L] [L'l ~ — / 'f^xipr) W[L][L']{r) jy{pr) dr + const , 
p J 
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and if the arbitrary constant of integration is taken to be we obtain for the short 

range behavior, 

1 f 

P[L\[u]i^^P) ,%^[L]W] - -] ^A(pr) WiL]iL']{r) J\'{pr)dr , (14) 

which obeys the condition (plsl). In practical calculations the last indefinite integral can be 
found analytically by using the leading terms of series expansions of nx, and jy. In 

this way one finds that the second term of Eq. (|14|) is regular at r = for A < A' (right 
upper corner and the diagonal of the matrix) and may be singular for A > A' (left lower 
corner of the matrix). However, this singularity is always compensated by the presence of 
jx in Eq. (|l^). Thus the coupled equations, Eqs. (||), along with the boundary conditions, 
Eqs. ([T3| ) and (p!^, form a well-defined differential problem. 

The regular basis obtained in the form of Eq. (^ is ideally suited for constructing 
physical solutions for one-channel problems. Indeed, since the right hand sides of Eqs. (^ 
vanish together with the potential, so do the derivatives drF^^\ which in turn implies that 
beyond some Vmax both function F\^^l'] l'^' P) become practically constant and the asymptotic 
behavior of is totally determined by the Riccati-Hankel functions. On the other hand, in 
one-channel problems we can have only three types of physical solutions describing bound, 
Siegert, and scattering states. For all of them the boundary function U of Eq. (H) must 
be constructed from the Riccati-Hankel functions which depend on the hyperradius and at 



large r behave as [12 



h^;^\pr) — i> =Fiexp[±i(pr )] . (15) 



Thus the physical solution is a linear combination of the form 



U[L]ir,p) = J2(l)[L][L']ir,p)A[L']ip) 

[L'] 



which at r = Vmax smoothly matches, and for r > Vmax automatically coincides with, the 
boundary function U, 



1 
2 



h^x\pr)Y. ^lL]lL']ir"^ax,p)A[L']{p) + h[ \pr)Y,F{L]lL']{rmax,p)A[L']{p) > = U[L]{r,p) 
\L'] \L'- 



provided the correct coefficients ^[l](p) are found. 

In the present work we are concerned with bound and Siegert states. For both of them 
each element of the column U[l] at large r must be proportional to h^^\pr) which expo- 
nentially decays when p is on the positive imaginary axis (bound state), or represents pure 
outgoing waves when p is in the fourth quadrant of the complex p-plane (resonant Siegert 
state). This can be achieved provided that 

T.^[L][L']irrnax,p)AlL']{p)=0. (16) 
[L'] 
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Here the argument of A is p because the wave function of a bound or Siegert state does not 
depend on the orientation of the incident momenta {p}. 



The homogeneous system of equations ([T6|) has a nontrivial solution if and only if 



detl|F[tj?^,j(w,p)l|=0. (17) 

The discrete points p = poi, po2, ... at which Eq. (|1^) is fulfilled, are the spectral points 
corresponding to bound and resonant states. 

When the number of particles N = 2, the functions F'^^\r^p) are closely related to the 
Jost solutions and the limit of F'^~\ 

f{p) = \im^F^-\r,p) , (18) 

is the Jost function. By analogy, we can call ||-Fj^^jj''^,j(oo,p)|| the Jost matrix for a one- 
channel A^-body problem. In practical calculations instead of r = cx) we can always use 



III. MATRIX ELEMENTS OF THE POTENTIAL 

/^From the simplest shell-model it follows that the lowest configurations for three and 
four neutrons are (Os)^(Op) and (Os)^(Op)^ and so we expect to find one and two neutrons 
with i = 1 respectively. More elaborate few-body analyses P,p!5|-[l9| corroborate this sim- 



plistic argument and concluded that if the nuclei and were to exist, the most favorable 
quantum numbers ( J'^) would be | for and 0+, 1+, or 2+ for ^n. This means that even if 
and do not exist, the resonant poles corresponding to these states must be the closest 
to the origin of the complex energy plane. In the present work we shall search only for these 
states. 

Since in the nnn and nnnn systems one and two particles respectively have ^ = 1, the 
minimal value of the grand orbital number for them is Lmin = 1 foi' and Lmin = 2 



for nnnn. The general rule defining Lmin can be found in Ref. |Tl|. It has been pointed 
out in many papers (see, for example, Refs. P, |l5| , ^ ) that the minimal approximation 
L = Lmin where only the first equation of the system (||) is retained, provides an adequate 
description of the properties of '^He and '^He nuclei and that the corresponding minimal 
components, U[L^.^]{r,p), of the wave functions contribute ~ 95% to the total normalization 
integral. Therefore we employ this approximation to investigate the analytical properties of 
the multi-neutron Jost function and thus to shed some light on the question of existence of 
resonances in many-neutron systems. 

The general form of the matrix elements W\t irr i for three and four neutron systems 
were given in Refs. 0] and |T^. They are 



<::'] (0 = f l^'^ de sin^ e cos^ e VUV2r COS 6) , (19) 
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(nnnn,l+) 

ra%n\ 



(n„nn2+) ,^ 



105-33 r/2 



8-16 Jo 

+^sin^^)Kn(V2rcos^) 



105-33 W2 



rf^ sin^ ^ cos^ B (4 cos^ ^ - 4 cos^ d sin^ ^ 



64 JO 



105-33 r/a 



sin^ d cos^ ^ 14n(v^r cos d) , 

Jo 



80 



/ rf^ sin^ 9 cos^ ^ (cos^ 9 + — sin^ ^)Kn( V2r cos 
JO 16 



These potentials are employed in our search for sub-threshold resonances. 



IV. COMPLEX ROTATION 



(20) 



(21) 

(22) 



Within the minimal approximation, only the first of Eqs. (0) remains and consequently 
we have only one pair of Eqs. (P). Therefore the problem is similar to the two-body one 
with potentials (p!9|-|22|) and angular momenta obtained by Eq. (^, that is A = 5/2 for {nnn) 
and A = 5 for {nnnn). The corresponding Jost function becomes an effective two-body one. 



It is known that starting with a two-body radial Schrodinger equation in its ordinary 
form, the Jost function for a long-range potential can be defined for Im {p} > only. The 
potentials (p!9|-p2D are clearly of the long-range type. This is a result of a general rule 
in which two-body potentials when sandwiched between hyperspherical harmonics, aquire 
slowly decaying tails (see Ref. [|11|]). For example, even if we take in Eq. (|1^) in the 
form of a square well, the resulting function M/^[L^i„][L^i„]('") at large r behaves as [16 



Thus, the limit (|T8|) does not exist in the fourth quadrant of the complex p-plane where we 
shall search for possible resonances. 

To overcome this difficulty we employ the complex rotation method in the form developed 
in Refs. In Eqs. (H) we replace the real hyperradius with a complex one, viz. 



a;exp(«t 



a; > 



< ^ < 



a,FW(x,^,p) 



d,F^-\x,e,p) 



2ip 



h'^^\pxe'^) F^+\x,9,p) + h[-\pxe''^) F^-\x,9,p) 



(-)/ 



^ •,{+) 



2ip 



hY\pxe'') W{xe 



h^^\pxe'^) F^+\x,9,p) + h['\pxe'^) F^-\x,9,p) 



(23) 



Such a rotation does not change the Jost function which is r-independent. Meanwhile, 
it changes the functions F*^^^ to the effect that F^~\r,p) can be defined above the line 
(— ooe~*^, +ooe~*^) in the complex p-plane. Moreover, in Ref. it was shown that at all 
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points above this line the hmit (0) exists and gives the correct Jost function. Therefore, 
using a rotation with large enough 9, we can calculate the Jost function at the points of 
interest in the fourth quadrant of the p-plane and thus the multineutron resonances, if any, 
can be located. 



V. RESULTS AND DISCUSSION 

A multitude of nn-interactions can be found in the literature. We may divide them into 
two classes, namely, those with and those without a repulsion at small distances. In this 
work we performed calculations using both types of potentials. 

As a purly atractive potential we choose the Yukawa type 

Vnnir) = -Vo , 

r 

with Vo = 54.7477 MeV fm and ao = 0.84034 fm^^. /^From the second class we employ the 
Malfliet-Tjon singlet potential ||2ll 



Vnnir) = V^—^-V2—^ , 

with Vi = 1438.72 MeV fm, V2 = 513.968 MeV fm, en = 3.11 fm"\ and as = 1.55 fm'^ 
Both these potentials describe a nucleon-nucleon interaction in the ^Sq channel and support 
a virtual di-neutron state. We note that we have slightly adjusted the depth of the Yukawa 
potential so that the virtual bound state it generates is the same as for the Malfliet-Tjon 
potential. This virtual bound state energy is given in Table 1. 

Solving Eqs. (^) numerically with the potential matrices (p!9| - p2|) , we found that, be- 
yond Xmax = 40 fm, the function F^~\x, 9,p) becomes practically constant (varying only in 
the 8-th digit). Further, the value of F^^\xmax,0,p) did not depend on the choice of the 
rotation angle 6 (indicative of the accuracy of the method) provided the point p is above 
the line (— ooe~*^, -|-ooe"*^). 

Using Newton's method, we have located one zero of F^~\xmax,(^,p) ioi each of the 
states ^n(l/2~), '^n{0'^), ^72(1"^), and ^n{2~^) in the fourth quadrant of the complex p-plane, 
for the Yukawa as well as Malfiiet-Tjon rm-potential. The coordinates Pq of these zeros are 
given in Table 1. All zeros found lie below the diagonal of the quadrant, which represents 
the threshold energy (ReE = 0); that is the corresponding energies, £"0 = pl/2m, have 
negative real parts. This implies that these zeros are sub-threshold resonances. 

It is evident from the results of Table 1 that the position of zeros are sensitive to the 
choice of the nn-potential. The movement of these zeros with an artificial increase of the 
potential strength also exhibits an extreme sensitivity to the choice of the potential. Thus, 
to generate and bound states, the Yukawa potential needs to be multiplied by 2.67 
and 2.32 respectively, while for the Malfiiet-Tjon potential by 4.04 and 3.59. 
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In Fig. 1 we show the movement of the three-neutron Jost function zeros on the mo- 
mentum plane, when the nn-potential is multiphed by a factor a, \4n — ^ CiVnn- Open and 
filled circles correspond to the Yukawa and the Malfiiet-Tjon potentials respectively. The 
zeros furthest from the origin correspond to the physical potentials {a = 1) and represent 
the zeros given in the Table 1. The upward sequences of points are shown for the uniform 
increase of a, a = 1.0, 1.1, 1.2, ■ ■ -. In Fig. 2 the corresponding trajectory in the energy 
plane is depicted. In Figs. 3 and 4 the results for the four-neutron system are plotted the 
notation being the same as for Figs. 1 and 2. It is worth mentioning that a trajectory similar 
to these figures, was obtained in Ref. for a sub-threshold resonance of the hypernuclear 
system Ann. 



There have been several previous theoretical attempts to locate three and four neutron 
resonances: using the hyperspherical approach p|, p!6|j23| -p5| , by an analytical continuation of 
the Faddeev kernel |]T^, and by using the complex scaling method Invariably the authors 
of these papers searched for multineutron resonances above the threshold and close to the 
real energy axis. Only one of these attempts turned out to be successful, the complex-scaling 
calculation by Csoto, Oberhummer, and Pichler where they found a three-neutron res- 
onance at E = (14 — ^6.5) MeV for the J'^ = I"*" state and that all other states up to J = | 
are nonresonant. 



In short, all the above authors agree that there are no resonances in the four neutron 
system, nor in the | state of three-neutrons. On the other hand, it has been shown many 
times that by an artificial increase of the neutron-neutron attraction, one can always obtain 
bound states for the ^n(| ) and ^?t,(0"'") systems. Since the corresponding poles cannot van- 
nish or undergo a discontinuous relocation, they must be smoothly move somewhere after 
decreasing the potential strength to its physical value. Thus, to say that the states ^n{^ ) 
and ^n(0"'") are nonresonant is unsatisfactory. One would like to know in which direction 
and how far the poles move. Our present work sheds some light on this problem. 

Using the minimal approximation L = Lminn means that the actual position of reso- 
nances in the complex momentum plane may be different from those shown in the figures. 
In this connection, however, we emphasize that in bound states calculations the minimal ap- 
proximation always underbinds the system and that the omitted higher harmonics amount 
to an additional attraction in the effective potential [^]. Since in this work we found that 
an increase in the attraction moves the resonances, in the energy plane, up (decreases their 
width) and to the right (closer to the threshold), the values given in Table 1 can be considered 
as lower bounds for the energies and upper bound for the widths of the resonances. 

Since the physical potentials generate negative real parts of i?o, the corresponding res- 
onances can be excited at negative energies only, that is when the or the system is 
placed in a sufficiently strong attractive external field. Such a situation can be realized, for 
example, inside a nucleus. Besides the negative energy, the interior of a nucleus provides 
a high frequency of multineutron collisions due to the high density of nucleons. Therefore, 
three or four neutrons inside a nucleus could form an unstable cluster corresponding to a 
sub-threshold resonance. 
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In conclusion, our approach which is based on the Jost function calculation, is a power 
method which enables us to investigate the analytical properties of the S-matrix in the 
complex momentum plane. This opens up new possibilities in locating sub-threshold reso- 
nances which is a difficult task for many other methods. Moreover the formalism given can 
be extended to complex values of the angular momentum A and therefore Regge trajectories 
can also be located. 
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TABLES 



potential 


Yukawa 


Malfliet-Tjon 




Po (fm ^) 


£^0 (MeV) 


Po (fm ^) 


Eo (MeV) 


nn (0+) 
nnn ) 
nnnn (0"*") 
nnnn (1"*") 
nnnn (2+) 


-i0.0401 
0.232 - i0.560 
0.337 - zO.594 
0.262 - i0.790 
0.467 - iO.744 


-0.0667 
-5.39 - i5.39 

-4.96 - i8.30 
-11.5 -i8.59 
-6.97 - il4.4 


-i0.0401 
0.291 - iO.569 

0.402 - iO.537 
0.435 - iO.783 
0.535 - iO.699 


-0.0667 
-4.95 - i6.87 
-2.64 - i8.95 
-8.77-il4.1 
-4.19-il5.5 



TABLE I. Jost function zeros in momentum and energy planes. 
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FIGURES 




FIG. 1. Movement of the 3n Jost function zeros when the nn-potential is multiphed by the 
enhancing factor Vnn — > oiVnn- Open and filled circles correspond to Yukawa and Malfliet-Tjon 
potentials respectively. The sequence of points are shown for a uniform increase of a by 0.1, 
a = l.O(lowestpoint), 1.1, 1.2, .... 



14 



Im E (MeV) 



-6 -5 -4 -3 -2 -1 

I I I I I L 



1 Re ^ (MeV) 
J I 



-1 - 



-2 - 



-3 - 



V 
: ° 

• o 



o 



o 



• o 

o 



o 



o -O 
o o o • 



-6 

-7 



FIG. 2. Movement of the three-neutron Jost function zeros on the energy plane. The notation 
is the same as in Fig. 1. 
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FIG. 3. Movement of the four-neutron Jost function zeros on the momentum plane. The 
notation is the same as in Fig. 1. 
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FIG. 4. Movement of the four-neutron Jost function zeros on the energy plane. The notation 
is the same as in Fig. 1. 
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